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Abstract 

We modify the Zee mass matrix by adding a real one parameter perturbation 
which is purely diagonal and trace-less. We show that in this way we can explain 
both solar and atmospheric neutrino oscillation data. There is a correlation between 
the deviation from strict maximality of |f/^3| = l/\/2, with the emergence of a 
small but non-zero C/eS- We calculate how big a value can [/es get when we restrict 
ourselves within the allowed regions of solar and atmospheric neutrino masses and 
mixing angles. We also discuss the impact of a 5*2 permutation symmetry on our 
mass matrix and show how a small t/es 7^ can emerge when this 5*2 permutation 
symmetry between the second and the third generation is broken. 
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1 Introduction 



With neutrino oscillations having been successfully established beyond all doubts by a series 
of spectacular experimental results |l|i2|i^i4j, the focus has now shifted to phenomenolog- 
ically determining the neutrino mass matrix and deciphering the underlying theory which 
gives us the correct neutrino mass matrix. Our current knowledge of the range of allowed 



values for the oscillation parameters at the 3a limit is [S] 

7.2 X IQ-^eV^ < Aml^ < 9.2 x IQ-^eV^ (1) 

0.25 < sin^ 9i2 < 0.39 (2) 

1.4 X lO^W < Aml^ < 3.3 x 10~W (3) 

sin^ 2^23 < 0.9 (4) 

sin^ 9i3 < 0.044 (5) 



A series of next generation oscillation experiments have been proposed/planned to mea- 
sure very precisely the oscillation parameters. Our knowledge on Am|]^(= Am^^^) and 
sin^^23(= sin^^ai^) is expected to improve to about 4.5% and 20% accuracy respectively 
[Hj, while that on Am2i(= AuIq) and sin^ ^i2(= sin^ 9q) could improve to 7% and 16% 
The precision on the mixing angle sin^ 9i2 could be improved further with new reactor ex- 
periments [ZllHl- The hitherto unknown mixing angle 6'i3 will be probed in the forthcoming 
long baseline and reactor experiments ^ to values of sin^ 6*13 ^ 10^'^ 6J. Any hint of CP 
violation in the lepton sector will be looked for in the future long baseline experiments. The 
deviation of ^23 from maximality and the sign of D23 = 0.5 — sin^ 6'23 can be experimentally 
checked in atmospheric neutrino experiments ^Ul- The sign of Arn^^ can be probed in 
either the long baseline experiments or in atmospheric neutrino experiments |TI] . It might 
even be possible to check the sign of Amg^ in experiments looking for neutrino-less double 
beta decay (Oz//?/5) *T^. The next generation OujSjS experiments are of course expected to 
improve the limits on the effective mass parameter to \{m) \ < 0.03 eV jTH]- Limits on the 
absolute neutrino mass coming from direct lab measurement is also expected to improve 
from its current limit of m^j < 2.2 eV to m/3 < 0.2 eV ^1]. The best current limit that we 
have on the absolute mass scale comes from cosmological measurements and these will be 
further improved in the future [T^\. Therefore we expect that the neutrino mass matrix 
will be determined fairly well by the next generation of experiments. 

In this paper we will phenomenologically analyze a modified version of the Zee mass 
matrix which is compatible with the current allowed values of the oscillation parameters. 
We will study the predictions this model makes for the oscillation parameters and check 
how it could be tested in the future. 
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The Zee ansatz is a theoretically attractive model of neutrino masses and mixings ^E] • 
This model needs a very minimal extension of the Higgs content beyond the standard model 
and one does not necessarily need to invoke supersymmetry. While the most general mass 
matrix in the Zee model has many more degree of freedoms than can be constrained by 
oscillation data, a more restricted texture for the Zee mass matrix emerges if one imposes 
a condition that only one of the two Higgs doublets in the model couples to the charged 
leptons J7j. This brings an added advantage that the resultant model does not suffer from 
problems concerning flavor changing neutral currents in the charged lepton sector. This 
ansatz leads to a neutrino mass mass matrix which is symmetric and for which the diagonal 
elements vanish. Therefore this mass matrix has only three real parameters. Consequently, 
when tested against experimental results, this ansatz falls in a situation which is extremely 
constrained [IHl EHl 120] • This is because one has to predict two mass differences and three 
mixing angles adjusting only three input parameters. It turns out that though Zee mass 
matrix can reproduce a maximal z/^ ^ Vr mixing very naturally, it fails to simultaneously 
reproduce the LMA region of solar neutrino oscillation. 

There has been many attempts to modify the Zee mass model [21]. There are mainly 
three points which are addressed in the literature which were looked into while Zee model 
was modified, (i) The Yukawa couplings which lead to masses of ordinary quarks and 
leptons are independent of the Yukawa couplings which generate the neutrino masses ra- 
diatively in Zee mechanism. There were attempts to link these two families of Yukawa 
couplings, (ii) There were attempts to embed the Zee model in grand unified scenarios, 
(iii) Zee model leads to severe constraint on sin^ 2^12 > 0.99. There were attempts to 
modify Zee model which will lead to correct values of sin^ 2^12 ~ 0.85 which is consistent 
with the LMA type solution. 

We modify the Zee mass matrix by adding a perturbation which is purely diagonal 
and trace-less. Furthermore the perturbation has only one real parameter. We do not 
introduce any specific field theoretic model for this extra diagonal elements. This extra 
piece is introduced in a purely phenomenological fashion. However, the sum of mass 
eigenvalues will remain zero as the mass matrix is trace-less. It has been shown from very 
general conditions that a trace-less mass matrix can fit the observed neutrino oscillation 
data well [22]- 

Therefore we consider a real symmetric trace-less neutrino mass matrix with four pa- 
rameters. We show that this real symmetric and trace-less four parameter mass matrix 
can correctly predict the mass squared differences and mixing angles needed to explain the 
world neutrino data. Because there are four input parameters and six testable observables 
are returned, one effectively ends up with two predictions. 

We begin in section 2 by briefly reviewing the problems faced and the current status of 
the original Zee-Wolfenstein ansatz for the mass matrix. We next analyze phenomenologi- 
cally the mass texture we obtain by modifying the Zee mass matrix by a trace-less diagonal 
perturbation matrix. We show that this mass matrix can return values of neutrino oscilla- 
tion parameters consistent with the world neutrino data. In section 3 we briefly discuss the 
models which could give our neutrino mass texture. We end in section 4 with conclusions. 
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2 The ansatz and the phenomenology 



2.1 The original Zee-Wolfenstein ansatz 



In the Zee model jTHI, under the assumption that only one of the Higgs doublets couple to 
the charged fermions ^^I, the neutrino mass matrix assumes the form 



(6) 



which is generally referred to as the Zee mass matrix in the literature. It can be shown that 
the three non-zero entries involved can be taken as real without any loss of generality and 
thus the Zee mass matrix is described by only 3 real parameters. Redefining sin^ = 
cos^ = and e = ^ we get 

mo ' mo 





( 






Mzee = 















/ 



/ sin^ cos 61 

sin^ e I . (7) 
Vcos^ e 



Mzee = mo 



Note that one of the main characteristics of this mass matrix is its trace-less condition 
Since the Zee matrix is real, this implies that the sum of its eigenvalues vanishes, 

mi + 1712 + "^3 = (8) 

This along with the condition that ee element of Mzee is exactly zero, leads one to conclude 
that the only neutrino mass spectrum allowed in this scheme is the inverted hierarchy JH] • 
If in addition we impose a. L^. — — lepton family symmetry j2Sl 121]; the entry 
e goes to zero and the resultant matrix reproduces bimaximal neutrino mixing as long as 
6' = 7r/4^. It also predicts ^m^^ = 0. To get a non-zero value for Amli one has to break 
the Le — Ln — Lr lepton family symmetry slightly, such that e is still smaller than the 
entries cos^ and sin^, with the difference between the respective magnitudes determined 
by the extent to which Lf. — L^ — is broken. However, even though this small breaking 
of the symmetry explains the solar mass squared splitting, it fails to drive the solar mixing 
angle far away from maximal mixing as required by the data. In fact, up to first order in 
the small parameter e, the solar mass splitting and solar mixing angle is given by 



Am^i ~ 2m^esin2^, (9) 



sin^ ^12 = ^ - ^ sin 2^. (10) 



^In fact it necessarily produces a maximal solar mixing angle and Ue3 — 0, while the atmospheric 
mixing angle is free to be determined by the value of 9 which could take any possible value. 
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This implies that the deviation of sin^ 612 from maximahty D12 = 0.5 — sin^ O12 is given by 

1^12 = J -^i? sin 20, (11) 

where we have used Amg^^ ~ m^. Since from the world neutrino data we expect the ratio 
R = Am2i/Am|^ to lie between [0.022 — 0.066], the maximum deviation of sin^6'i2 from 
its maximal value that we can have in this model is 0.00825. This is in stark contradiction 
to the solar neutrino data, which predicts a deviation of at least D12 > 0.11 at the 3a CL. 
In fact, the condition D12 < 0.00825 is ruled out at more than 6a from the solar neutrino 
data - meaning that the simple Zee-Wolfenstein ansatz is disfavored by the data at more 
than 6a CL. 



2.2 The perturbed Zee ansatz 



We add to the mass matrix Eq. (jU)) another perturbation matrix which is diagonal and 
trace-less. 

/-2 m 



Redefining sinO = cos0 = e = ^ and 5 = ^ we get, 

mo ' mo ' mo mg ° ' 





( 




M = 
















'12) 



M = niQ 



—26 sin 6 cos 6' 
sin 6 6 e 
cos 9 e 6 



(13) 



An additional u,. 





1 


-26 


1 




M' = niQ 




1 


6 


6 




V 


1 

V2 


6 


^ J 



symmetry [23] will force 6 = n/A and simultaneously e = 6 giving 



(14) 



This form of the mass matrix has been independently derived from SU(3) global flavor 
symmetry j^Sj. However, in this case t/es is strictly zero and 623 is strictly maximal^. We 
will work in a scenario where e can be different from 6, while both are kept small compared 
to cos 6 and sin 6. We will also let 6 to take any possible value. This would results in 623 
deviating from maximal and Ues from zero and would enable us to relate the deviation of 
9i2 from its maximal value to explain the LMA solution, to the deviation of 623 and Ue3 
from their maximal and null values respectively. 

In the limit that 6 and e are considered small, we keep only up to the first order terms 
in these parameters and obtain ^ 



Am^^ ~ m^(-25 + 2esin20) 



(15) 



^We will discuss this case again towards the end of this section. 

■^Note that the neutrino mass hierarchy predicted in this model is inverted and sgn^Am'^i) is negative. 
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|Am3i| ^ ml , (16) 
sin'^ia ~ i - i(35 + esin2^) , (17) 

tan^^,3-^, (18) 

Ue3^V2ecos29. (19) 

From these expressions we note that the deviation of 623 from its maximal value is directly 
related to the deviation of Ues from zero and in this model the two are related through 

Requiring that 0.022 < R < 0.066 {R = Am2i/ Am"^^) to satisfy the experimental data, we 
derive the condition using Eqs. (fT3j) and (fTBj) 

0.011 < (-5 + e sin 2^) < 0.033 . (21) 

We further note that neither 6^23 nor Ues are dependent on the value of 6 and depend on 
e and 6 only. The solar mixing angle on the other hand is given from a combination of 6 
and e. From the condition that sin^ 612 < 0.39 at 3cr, we derive that 

35 + esin2e > 0.44 . (22) 

One can check that Eqs. (|2H) and (j22j) are completely compatible. This implies that our 
ansatz for the neutrino mass matrix should be completely compatible with world neutrino 
data. Combining Eq. (j^^ with Eq. (PT|) we get the approximate conditions 

0.107 <5 < 0.242 , (23) 
0.119 0.275 

We show in Fig. |21the values of the parameters e and S for which our mass matrix predicts 
values for the oscillation parameters consistent within their current 3a allowed limits (given 
in Eq. (Q)-©). The parameter 6 of the mass matrix is allowed to vary freely and we show 
the results for three different values of mo given in the caption of Fig. |21 The allowed 
values for 6 and e roughly are seen to correspond to those given by Eqs. (j23|) and (I24j) . In 
Table 1 we give the values of the oscillation parameters predicted by our mass matrix for 
a fixed set of mo, 6, e and 6. 
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If we denote the deviation of 9^ from its maximal value by D12 = 0.5 — sin^ 6*12 and 
deviation of 923 from maximality by D23 = 0.5 — sin^ 923, then 

1 



D,2 = -[36 + e^l-ADl,j , (25) 

while the deviation of Ues from the value of zero is given by 

U,3 = 2V2eD23 . (26) 

We show in Fig. El the ranges of oscillation parameters predicted by our neutrino mass 
matrix. In this figure we fix ml at 0.0015 eV^ (red dots), 0.002 eV^ (green dots) and 0.003 
eV^ (blue dots) and let e, 6 and 9 take on any possible value such that the world neutrino 
data is satisfied at the 3a CL. The left panels show the correlation between the predicted 
mixing angles while the right panels give the variation of the ratio R with the mixing 
angles. Note that while the predicted mixing angles and R themselves are independent of 
rriQ, the apparent dependence seen in the figure comes from the fact that the mass matrix 
has to simultaneously explain the individual mass squared differences Am^i and Am|]^, 
which depend on mg. 

We next turn our attention to the predicted value for the effective mass | (m) | that can 
be observed in neutrino-less double beta decay (Oz//?/?) experiments. Since this is given 
by just the modulus of the ee element of the mass matrix, the predicted effective mass in 
Oiy/3/3 is 

\{m)\=2mo5. (27) 
The prediction for square of the mass observed in tritium beta decay experiments is 



771 

^__0[2 + esin2^], (2J 



2 



while the sum of the masses which is constrained from cosmology is predicted in our model 
to be 



\mi\ ~ 2mc 



1 + ^(6 -esm29) 



(29) 



The three observables which depend on the absolute value of the neutrino masses are 
related in our model as 

, , I ("m) I 2ml , ^ 

~ 2 mo 

Note that while \{m)\ depends on 6, m| on e and 9 and Y.i on e, S and 9, their 
combination Eq. (jHU)) depends only on mo. We show in Fig. El the predicted values for 
the observables in tritium beta decay (m/3), Ou^P (K"^)|) and cosmology (X^i |"^j|) in the 
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Texture parameters 


sin^ 


2^0 


Sin^ 29atm 






Ue3 




mo 


e 


6 


e 








(eV^) 


(eV^) 




1 


0.039 


0.21 


0.21 


7r/4 


0. 


85 


1.0 





1.8 X 10~3 





2 


0.039 


0.21 


0.236 


0.84 


0. 


82 


0.99 


9.0 X 10-^ 


1.8 X 10-3 


0.019 


3 


0.039 


0.21 


0.236 


0.73 


0. 


82 


0.99 


9.0 X 10^^ 


1.8 X 10-3 


0.019 



Table 1: Some specific values of input parameters and outputs. In figures more details can 
be found. It can be noted that atmospheric mixing angle deviates from strict maximality 
when a non-zero Ug^ is returned. 



left panels. The right panels show the correlation between the effective mass | (m) | and the 
oscillation parameters. 

We reiterate that our mass matrix has four input parameters namely mg, S, e and 
6. The overall mass scale as well as the mass squared differences are controlled by mo, 
whereas it does not influence the mixing angles at all. Thus mo, which is constrained by 
the solar and atmospheric mass squared differences will have an upper bound coming from 
experiments on beta decay, neutrinoless double beta decay and cosmology. The mixing 
angles are controlled by the three other parameters, namely, S, e and 6. If we impose the 
condition of maximal 623 on the mass matrix, then we must have 6 = tt/A (cf. Eq. ()18|)). If 
we had a S2 exchange symmetry between the second and the third generation (z/^ ^ z/^), 
we will get 6 = tt/4. However, this symmetry also makes e = 6 and we get the mass matrix 
given by Eq. (jl4|) . The predicted values of the oscillation parameters for this case is given 
by the flrst row of Table 1. As discussed before, in this case the mixing angle 623 is exactly 
maximal, while Ue3 is exactly zero. However, note that the solar mass squared difference 
Amgi is also exactly zero. That the predicted 623 is maximal, Ues = and Amgi = in 
this case, can also be seen from Eqs. p8|) . (fT^ and (fT3j) respectively. This mass matrix 
is therefore phenomenologically untenable. In order to be able to explain the neutrino 
oscillation data, this exchange symmetry will have to be broken. This discussed in 

detail above, is consistent with all observations. The braking of the symmetry forces S ^ e 
giving a small Am^i and 6 ^ n/4 causing 623 to deviate from maximality and Ue3 from 
zero - all of which are naturally small, being protected by the approximate 5*2 symmetry. 



3 Discussions on possible models 

Let us write Eq. (fT^ as 

M^j = Xij + Yij (31) 

Here Xij is the symmetric off-diagonal part and Yij is purely diagonal and trace-less. In 
this paper we work with SU{?))f symmetry as was used in |26j but we consider a further 
step. We introduce SU{?))f symmetry breaking also; which is obtained by the VEV of an 
octet Higgs ifg 

SU{'i)F ^^^^ SU{2)f X U{1)f (32) 
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Leptons are triplets and antileptons are antitriplets whereas quarks are singlets. When 
SU{3)f symmetry is broken, electron generation becomes a singlet whereas second and 
third generation remains as a doublet of SU{2)p. This helps maximal mixing of atmo- 
spheric neutrinos. 



(33) 
(34) 



It is easy to see that this symmetry breaking pattern will be obtained if Hg gets a VEV of 
the form 




< H. >-- 




(35) 



For explaining this let us write down how SU (3) p representations transform under SU{2) p x 
U{l)p. 



3 (1,-2) + (2,1) 

3 (1,2) + (2,-1) 

8 (2, 3) + (2, -3) + (1,0) + (3,0) 



(36) 
(37) 
(38) 



Neither 3 nor 3 can get VEVs because then we will also break the residual symmetry 
SU (2) p X U{1) p. The (1, 0) component of 8 must get the VEV; because that will leave the 
residual SU{2)p x U{l)p intact. In matrix form the (1,0) component of the octet Hg is 
written in Eq. EHl The 2x2 submatrix in the lower right hand corner is a unit matrix so 
this 2x2 submatrix is a SU{2) singlet. The trace is vanishing so U{l)p quantum number 
is zero. Thus this specific symmetry breaking pattern fixes the form of VEV given in Eq. 

ESI 

From the form of the VEV we can see that a residual SU{2)p symmetry remains among 
the second and the third generation which assures maximality of mixing. The list of leptons 
Higgs scalars and their SU{3)p x SU{3)c x SU{2)p x U{1)y representations are 

FERMIONS 



L 



(3,1,2,-1/2) 
(3,1,1,1) 



(39) 
(40) 



ZEE MODEL SCALARS 
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00 = (1,1, 2, -1/2) (41) 

01 = (3, 1,2, -1/2) (42) 
X= (3, 1,1,1) (43) 
77 = (1,1, 2, -1/2) (44) 

EXTRA SCALARS 

Al = (3,1,3,1) (45) 

i^8 = (8,l,l,0) (46) 

(47) 



First let us explain Yij term of Eq. (jSH). The Higgs field is a SU{2)l x U{1)y 
triplet as well as a fiavor triplet which couples to z/j and uj and gives a direct left handed 
Majorana mass term of the form fijViVjVLi where vl is very small. This mass term acts 
as the diagonal and trace-less perturbation if the Yukawa coupling fij is diagonal and 
trace-less. 

Let us explain why the coupling fij is of this specific form which is diagonal and trace- 
less. The reason is that it has its origin in a higher dimensional operator of the SU (3) 
fiavor group, 

^(3x3x3x8) = j^LLAl < Hs > (48) 

At this stage the coupling / does not have any generation index. Putting the value of 
< Hs > from Eqn we get the generation indices as, 

h.UL^AL (49) 

Now fij retains the symmetry between the second and third generation as well as it is 
diagonal and trace-less. 

In the Xij term of Eq. ()31|) is symmetric with diagonal entries vanishing. It can be 
generated by the following diagram in Fig ^ in a Zee type model in the presence of our 
SU{'i)F fiavor symmetry. The VEV of 0i(3) is of the order of electroweak scale. This 
diagram has been studied very well in the literature. The off diagonal mass matrix is of 
the form 

niij oc gij [m] - ml) (50) 

Here gij is an off-diagonal Yukawa coupling El • Actually gij controls the coupling 
strength between two lepton doublets and the charged Higgs Fig H Therefore 

SU{2)l antisymmetry forces it to be off-diagonal. So, when we multiply gij by (m^ — mf) 
the product becomes symmetric with vanishing diagonal elements. 
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(8) (01 (3)) 



▲ 



/ \ 
/ \ 



> • <- — X ► » 



T 



Figure 1: Extra SU{2)i doublet 0i transform as 3 of SU{?>) flavor. It does not couple 
to fermions at tree level. Quarks do not feel SU(3) flavor in this model. Leptons and 
antileptons are triplets and antitriplets of SU(3) flavor, x is a SU{2)l singlet whereas rj is 
a doublet as in minimal Zee model. In bracket we have displayed the SU{3) flavor quantum 
numbers. 

Many different models, in principle, can be constructed which are aesthetically more 
appealing than this, which will lead to Eq. (jl2j) . We do not know which one is right 
and which is wrong. Typically the first part may be generated by a modified Zee type 
mechanism and the second part may be generated by a Higgs mechanism or a see-saw 
mechanism or vice- versa. In this paper we do not focus on the details of the model building 
apart from citing an example. More models and details based on various flavor groups, 
will be presented in a future paper [27]. Here we have added the diagonal and trace- less 
perturbation on a purely phenomenological basis. 

4 Conclusions 

In this paper we have tried to generate simultaneously a non-zero t/es and LMA solar 
neutrino mixing angle in a Zee type model by putting in a diagonal and trace-less per- 
turbation. While doing so we have succeeded to keep Am^i and Am|^ as well as the 
atmospheric mixing angle within experimentally allowed ranges. The resulting mass ma- 
trix texture is symmetric and trace-less. It has four input parameters and six outputs 
therefore giving two predictions. When atmospheric neutrino mixing angle is strictly max- 
imal, [/e3 vanishes. However when atmospheric neutrino mixing angle deviates from strict 
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maximality, a non-vanishing Ues emerges. We studied the correlations between the dif- 
ferent oscillation parameters as well as observables depending on the absolute neutrino 
mass scale, such as the effective mass in neutrino-less double beta decay, mass parameter 
relevant in beta decay and the sum of the neutrino masses relevant for cosmology. The 
neutrino mass hierarchy is predicted to be inverted and the oscillation parameters given 
by this mass matrix are well within the reach of future experiments. 
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Figure 2: The values of e and 6 which satisfy the current 3a hmits from the world neutrino 
data. The red dots correspond to mg = 0.0015 eV^, the green dots to ml = 0.002 eV^ and 
the blue dots to mp = 0.003 eV^. The parameter 6 of the mass matrix is allowed to take 
any possible value. 
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Figure 3: Scatter plot showing the correlation between the values of the oscillation param- 
eters allowed by our perturbed Zee mass matrix. The red dots correspond to ml = 0.0015 
eV^, the green dots to ml = 0.002 eV^ and the blue dots to ml = 0.003 eV^. 
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Figure 4: Left columns give the scatter plots showing the correlation between the values of 
< m >, and allowed by our perturbed Zee mass matrix. Right columns show the 
correlations between < m > and the mixing angles and ratio R. The red dots correspond 
to ttIq = 0.0015 eV^, the green dots to ml 
eV^. 
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